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Abstract 

In this paper we report the experimental study of high-pressure xenon used 
as a scintillator, in the context of developing a gamma ray detector. We 
measure a light yield near 2 photoelectrons per keV for xenon at 40 bar. 
Together with the light yield, we also measured an energy resolution of ~9% 
(FWHM) at 662 keV, dominated by the statistical fluctuations in the number 
of photoelectrons. 
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1. Introduction 

Xenon detectors have been used and are in use for a number of applica- 
tions, in particular for the detection of gamma rays (for a partial list see [1] 
and references therein, For exhaustive review work see P| and [7]). 

Xenon detectors are also preeminent in the field of Dark Matter searches and 
neutrinoless double beta decay (e.g. the XENON experiment [8j, the LUX 
experiment [9] , and the EXO experiment [10] ) . Xenon is an attractive mate- 
rial for gamma ray detection, in particular thanks to its high atomic number 
(Z = 54) and, as a consequence, large cross-section for photoelectric absorp- 
tion. Xenon is also a rather dense material (the density of liquid xenon is 
3.1 g/mL, while the density of 50 bar xenon at 293 K is 0.4 g/mL), which 
allows to build compact and efficient detectors. The scintillation light from 
Xe is at 175 nm, in the VUV, and it is therefore rather difficult to detect 
efficiently. While a number of works presenting high-pressure Xe ionization 
chambers can be found in the literature, very few are available that discuss 
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the properties of high-pressure Xe as a scintillator (known to these authors: 
[11]). Compared to a liquid Xe detector, a high-pressure Xe detector does 
not require cryogenics, and this simplification may be important for some 
applications. It has to be noted that the energy resolution achieved in Xe 
gas proportional chambers is significantly better than the one in liquid Xe 
|12j . It is therefore very interesting to compare the energy resolution ob- 
tained from primary scintillation in Xe gas to the one measured for liquid 
Xe. 

In this work we have performed an experimental study of the properties 
of high-pressure Xe as a scintillator, in the context of developing a gamma 
ray detector for the detection of Special Nuclear Material^ Our first goal 
has been to study experimentally the light yield and energy resolution from 
a high-pressure Xe tube, studying the dependence on thermodynamic condi- 
tions as well. A short description of the apparatus is given in Sec. [2| and the 
results are discussed in Sec. [3l 



2. Experimental apparatus 

We have used a high-pressure tube (designed to withstand pressures up to 
200 bar) filled with high purity xenon [13j. The tube has a length of 200 mm 
and a diameter of 44 mm, with UV transparent windows at both ends. The 
inner surface of the tube is lined with a refiector combined with a wavelength 
shifter, that shifts the VUV light from xenon scintillation into visible lighij^ 

The isochoric curves of xenon are shown in Fig. [T] The critical point for 
xenon is at 289.77 K and 58.41 bar, therefore at room temperature (293 K), 
xenon is in supercritical phase (for pressures above 58.41 bar), while for 
temperatures below 290 K it becomes liquid. The detector was operated in 
the pressure range between 30 and 60 bar, i.e. with xenon as a dense vapor 
(density changing between - roughly - 0.2 and 0.7 g/mL). 

The scintillation light from the high-pressure tube was read out using 
two Hamamatsu R580 (quartz window) photomultiplier tubes (PMTs). A 
schematic diagram of the experimental setup is shown in Fig. |2} The data 
were taken using a CAEN V17510 digitizer (1 GS/s, 10 bit). A waveform 



^Project MODES_SNM, IhttpV/www.modes-snm.eu/' 

^The tube is produced by Arktis Radiation Detector Ltd. (http://www.arktis- 



detectors.com/[ ) 

''iittp://www.caen.it/csite/CaenProd.jsp?parent=ll&idmod=602 
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Figure 1: Isochoric curves for xenon. Data from NIST 
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is shown in Fig. [3j The two PMTs were digitized independently, and the 
trigger, directly programmed on the digitizer, required the coincidence of the 
two PMTs. 

We took data with various gamma ray sources, most of the time with the 
source coUimated with lead bricks, having a 3 mm wide slit illuminating the 
detector. 

2.1. Analysis procedure 

Our analysis procedure consists of four steps: 

1. Single photoelectron (p.e.) calibration: The single p.e. is taken from 
the digitized events. A peak finding algorithm is applied to the digi- 
tized waveform (Fig. [3]) and the area for each peak (in units of ADC 
counts xns) is determined integrating the signal over the time window 
in which the waveform is above threshold; the threshold is 2 times the 
noise (rms of the baseline), and a minimum of 3 consecutive samples 
above threshold is required. The time windows for each peak are shown 
in Fig. [3j The interesting region for the single p.e. is easily determined, 
producing a spectrum as the one shown in Fig. |4j The superimposed 
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Figure 2: Schematic diagram of the experimental setup. See text for explanations. 



function, fitted to the spectrum, is 

n=3 

Aog{0, s) + ^ Ang{n x m,y/n x a), 

n=l 

where Aog{0, s) is a Gaussian with amphtude Aq, mean zero and sigma 
s to account for the background, and the sum of the Gaussian functions 
models the photoelectron peaks. 

2. Each waveform is analyzed: The procedure is similar to the one used 
for the single p.e. calibration, but this time the amplitudes of the peaks 
in the event are summed up, to define the number of p.e. per PMT for 
the given event. This integration procedure is performed over a time 
interval of 400 ns following the first threshold crossing. 

3. At this point it is possible to apply basic selections to the reconstructed 
events. We have three different selections: first of all, the time differ- 
ence between the first identified peaks of the two PMT signals must 
be less than 40 ns (this selection rejects less than 0.1% of the events). 
Second, the difference between the integrals of the two PMT signals is 
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Figure 3: Digitized waveform from one the PMTs. The samphng rate is 1 GS/s, with a 
resolution of 10 bit. Also shown are the integration windows for the signal (see text for 
explanation) . 



required to be less than 200 p.e.. Third, the integrals of the two PMT 
signals are required to differ less than 30%. The last two selections on 
the asymmetry between the two PMTs tend to improve the peak-to- 
total ratio in the energy spectrum and are especially effective for data 
with the source collimated in the center of the detector, because events 
showing largely asymmetric signals between the two PMTs tend to be 
only partially contained. One example showing the impact of these se- 
lections is given in Fig. [sj for a collimated ^^''Cs source: after applying 
the selections described here, the peak-to-total ratic|^goes from 3.8% 
to 7%. 

The event integrals for each PMT are summed up, giving the amount 
of light for each event (in units of p.e.) and the final spectrum. 



"'Here for the "total" we integrate over the entire spectrum, including the X-ray lines. 
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Figure 4: The single p.e. obtained from the digitized events. The superimposed function, 
fitted to the spectrum, includes a Gaussian to account for the background, and a sum of 
Gaussian functions to model the peaks of the n p.e. (n — 1, 2, ...). 

3. Results 

In the development of a gamma ray detector, our main goal was the 
determination of the light yield of high-pressure xenon, when exposed to 
gamma rays. 

The energy spectra from exposure to a ^^^Ba and a ^^^Cs are shown in 
Fig. [6] (left and right, respectively). In both cases the source was coUimated 
at the center of the detector by lead bricks. ^^^Cs gives a line at 662 keV 
(and 32 keV X-ray line), while ^^^Ba gives several lines: the most prominent 
are at 31 keV (X-ray), 80 keV, 302 keV and 356 keV, with additional weaker 
lines at 276 keV and 384 keV, which, with the available energy resolution, 
are partially merged with the 302 keV and 356 keV lines, and hard to detect 
individually. The backscattered peak is large for this particular detector 
geometry, where the mass of the stainless steel tube largely exceeds the mass 
of the xenon inside. Data were taken also with a ^^Na source, which gives 
lines at 511 and 1275 keV. From the data, we derive a light yield of 1.91 ± 
0.05 p.e./keV at 662 keV. 

The energy dependence of the light yield (per unit energy) is shown in 
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Figure 5: Energy spectrum of a coUimated '^'^''Cs before (continuous red line) and after 
(dotted black line) the selections described in the text. 

Fig. [7| where we present the "normahzed hght yield" , i.e. normahzed to the 
hght yield per keV at 662 keV (^^"^Cs). There is a systematic trend, showing a 
larger light yield at lower energies. We interpret this as due, at least in part, 
to non-containment of the electron from gamma conversion. Fig. |8] shows 
the fraction of energy deposited in a pressurized xenon tube as a function 
of electron energy, for different gas pressures. Here monochromatic electrons 
(1000 for each energy) are generated in the center of a tube with a diameter 
of 40 mm, using the GEANT4 simulation packag^ At low pressures and 
densities the fraction of energy deposited drops rapidly with energy, because 
MeV electrons are, on average, not contained. At higher pressures the elec- 
trons stop inside the tube, but energetic Bremsstrahlung photons still escape 
from the tube. Considering the CSDv^ range for electrons, we see that it 
exceeds 20 mm in 40 bar xenon when the energy of the electron is larger than 
1.1 MeV: given a diameter of the detector of 44 mm, the probability of full 
containment of electrons with a range exceeding 20 mm is rather small. 



^http: / / geant4.cern.ch/ 

^Continuous Slowing Down Approximation, which gives an upper limit to the range 
of the electron. It can be easily calculated using the ESTAR program from NIST, 
http://physics.nist.gov/PhysRefData/Star7Text7MTAR:html[ 
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Figure 6: Left: "3Ba spectrum. ^^^Ba gives lines at 31 keV (X-ray), 80 keV, 302 keV 
and 356 keV. At ^30 keV several X-ray lines from iron and lead are also possible. Right: 
^^^Cs spectrum, no off-line selections applied. Both sources were collimated in the center 
of the detector. 



Again in Fig. [7], the measured energy resolution and the expected energy 
resolution (for the measured light yield, and assuming that the spread can be 
modeled with a Poisson distribution), are shown. The measured energy res- 
olution is consistently worse than the one expected from simple photostatis- 
tics. Again, the discrepancy increases with energy. The non-containment of 
the electron is one reason for this systematic effect. Another reason is that 
Compton scattering becomes more and more important when the energy of 
the gamma ray increases, and therefore the event is less localized and the 
energy is deposited within a larger fraction of the active volume. This leads 
to a degradation of the energy resolution, as can be seen in Fig. [9] and in 



Fig. 10, which show that in our setup the light yield has a significant (of 
the order of 15% over the full volume) dependence on the interaction posi- 
tion inside the detector. In Fig. |9] we show the number of photoelectrons 
collected on one PMT vs. the one from the other PMT, for an uncoUimated 
^^^Cs source at a distance of 2 m. It is clear that the band corresponding 
to the full energy peak has a "banana shape", i.e. when the sharing of the 
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Figure 7: Normalized light yield, energy resolution, and limit to the energy resolution 
from photostatistics, vs gamma ray energy. The normalized light yield is defined as the 
ratio between the light yield at a given energy divided by the light yield at 662 keV (^^^Cs 
line). The energy resolution is derived from fitting the peaks in the spectrum with a 
Gaussian plus a second order polynomial for modeling the background. The "statistical 
resolution" is calculated from the number of photoelectrons, assuming that the spread in 
the peak can be modeled with a Poisson distribution. 



light is almost equal between the two PMTs, the total collected light is less 
than when the light distribution between the two PMTs is very asymmetric. 



The data shown in Fig. 10 were taken in different conditions, with a coUi- 
mated ^^^Cs source, scanning the length of the tube up to 7 cm from the 
tube center. The dependence of the number of photoelectrons with respect 
to the position of the source is quite clear. In fact, even when summing up 
both PMTs, we notice that the light yield for an event in the center of the 
detector is about 15% lower than for an event happening next to one of the 
two PMTs. To summarize, the lower the energy of the gamma ray, the better 
the localization, because 



1. the range of the secondary electron is shorter; 

2. the gamma ray is more likely to deposit all the energy in one single 
interaction (photoelectric absorption, to be contrasted with multiple 
Compton interactions). 

A better localization immediately translates in a reduced spread on the light 
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Figure 8: Fraction of energy deposited in a pressurized xenon cylinder of 20 mm radius for 
electrons starting at the center of the cylinder as a function of the electron initial energy, 
for five different pressures. 

yield, and on the energy resolutioiiQ 

With the help of a Monte Carlo simulation we estimate the light collec- 
tion efficiency of 17%, not including the quantum efficiency of the PMTs, 
with large systematic errors due to the lack of independent knowledge of the 
optical properties of the various materials. Given the geometry of the tube 
and an aspect ratio of almost 5, only a small fraction of the scintillation light 
is directly detected by the PMTs, without being converted after interaction 
on the detector walls. The light collection efficiency takes into account the 
propagation of the light in the tube, i.e. the wavelength shifting process on 
the detector walls (for VUV photons: assumed to be 100% as estimated from 
available data), the reflection from the detector walls (for visible photons), 
the transmission/reflection from the windows (separately for VUV and visi- 
ble photons), and the possibility of absorption. The poor knowledge of the 
reflectivity of the material lining the tube introduces the largest uncertainty 
in our calculation. If we keep the reflectivity as the only free parameter in 
the calculation, we can determine it from the data, in particular from the 



partial correction for the position dependence of the light yield can be introduced 
offline, allowing to recover the energy resolution. 
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position dependent light yield as shown in Fig. 10 We estimate a reflectivity 
of 93% ± 2%, which translates into a light collection efficiency of 17l|% The 
quantum efficiency of the PMTs is ~27%, both for visible and VUV pho- 
tons. Therefore we measure an average energy expenditure (Wg) of 26lg eV 
to produce one scintillation photon in high-pressure xenon at 662 keV. This 
can be compared to a value of 23.7 eV for liquid xenon |14j . 

For completeness, we notice that our value for Wg is significantly smaller 
than the one of 76 ± 12 eV from [15] for 15 bar Xe at 60 keV, and 72 ± 6 eV 
reported in [16] at much lower pressure, between 1 and 3 bar, at 5.9 keV. An 
even larger Wg of 111 ± 16 eV was reported in [T7j for 1 bar pressure at 5.9 
keV. When comparing these conflicting results one should also consider the 
different thermodynamic conditions of the xenon. 




collected light - left - (p.e.) 



Figure 9: PMT-left vs PMT-right. Detector exposed to an uncoUimated ^"^ Cs source. 



We have also studied the light yield as a function of the xenon pressure: 



results are shown in Fig. 11 From the point of view of an application, this 
is important because an efficient detector should be operated at the largest 
possible density. We clearly see that the energy resolution deteriorates when 
the xenon gas gets closer to the supercritical phase, essentially because many 
scattering centers within the fluid hamper the uniform propagation of light. 
The abrupt change in the optical properties of the Xe is easily observable by 
simply watching the Xe cell while changing the pressure: the high pressure Xe 
goes from a perfectly transparent medium to a turbid one when the pressure 
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Figure 10: Collected light (normalized) as a function of the source position. Detector 
exposed to a coUimated ^^^Cs source, moved along the detector. 



exceeds 45 bars. In fact, while the hght yield does not change enough to 
modify much the Poisson term in the resolution, the effect is particularly 
visible on the measured energy resolution, pointing to a systematic smearing 
which applies to the light production or to the light propagation, or to both. 

It is interesting to compare our results in terms of energy resolution (9% 
FWHM at 662 keV, with a light yield of 1.9 p.e./keV) to the results ob- 
tained using liquid Xe detectors: Akerib et al. |9] report a light yield of 8 
p.e./keV and an energy resolution of 11.7% FWHM at 662 keV, while Ni et 
al. [H] report 6 p.e./ keV and 14.5% (FWHM at 662 keV). Therefore we 
see that, while the energy resolution in high-pressure xenon is dominated by 
photostatistics, this is not the case for liquid xenon. 

4. Conclusion 

In this paper we have presented experimental results that give an initial 
characterization of the properties of high-pressure Xe g clS clS db scintillator, 
in particular in terms of the light yield. Using a detector which is not fully 
optimized for efficient light collection, we measure a light yield of 1.91 ± 
0.05 p.e./keV at 662 keV at 40 bar pressure. The corresponding energy 
resolution is 9% FWHM. We find this result encouraging for applications of 
high-pressure Xe as a scintillator for gamma ray spectroscopy. 
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Figure 11: Light yield and energy resolution as a function of the xenon pressure. Detector 
exposed to a coUimated ^^^Cs source. 

Some of the shortcomings of the present setup are well understood in 
terms of non-uniformity of the light yield and can be readily fixed in an 
optimized setup. A more uniform light collection across the active volume 
should push the energy resolution down to the limit given by photostatistics, 
which right now is reached only at the lowest energies. We have also studied 
the behavior of xenon as a function of the pressure, and we find that, when 
the gas gets close to its supercritical phase, the energy resolution degrades 
significantly, well beyond the effect of the slightly reduced light yield. This 
systematic effect is due to the presence of many scattering centers that affect 
the propagation of light in xenon, and determines an upper limit to the 
density of xenon used as a scintillator. 
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